Scalar field cosmology: its non-linear Schrodinger-type formulation 
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Scalar field cosmology is a model for dark energy and inflation. It has been recently found 
that the standard Friedmann formulation of the scalar field cosmology can be expressed in a non- 
linear Schrodinger-type equation. The new mathematical formulation is hence called non-linear 
Schrodinger (NLS) formulation which is suitable for a FRLW cosmological system with non-negligible 
barotropic fluid density. Its major features are reviewed here. 
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I. INTRODUCTION 



The present universe is under accelerating expansion. This is convinced by many present observational data from 
cosmic microwave background Q , large scale structure surveys [3] and supernovae type la [1, [3, @ ■ There are many 
— i' ideas to explain such an expanding state, mainly it can be classified into three types: braneworlds and modification 
^ I of gravitational theory (e.g. Q), backreaction effect from inhomogeneity 0] and dark energy (for review, see [8|]). 
I ^ Dark energy is a type of cosmological fluid appearing in the matter term of the Einstein equation with equation of 
^ . state WD.E. < —1/3 so that it can generate repulsive gravity and therefore accelerating the universe. The simplest 
dark energy model is just a cosmological constant with = — 1. However the cosmological constant suffers from 
fine-tuning problem. Observational data suggests that the present value of wd.e. is very close to -1 and it also allows 
possibility that dark energy could be dynamical. Therefore scalar field model of dark energy became interesting 
^ ' topic in cosmology since time-evolving behavior of the scalar field gives hope for resolving the fine-tuning problem. 
Although the scalar field has not yet been observed, it is motivated from many ideas in high energy physics and 
quantum gravities. Theoretical predictions of its existence at TeV scale could be tested at LHC and Tevatron in 
. very near future. Phenomenologically the scalar field is also motivated in model building of inflation where super-fast 
acceleration happens in the early universe 0] . Cosmic microwave background data combined with other results allows 
possibility that the scalar field could be phantom, i.e. having equation of state coefficient < — 1 [lo| . The phantom 
equation of state is attained from negative kinetic energy term in its Lagrangian density The most recent 

', five-year WMAP result [l^ combined with Baryon Acoustic Oscillation of large scale structure survey from SDSS 
and 2dFGRS fli| and type la supernovae data from HST 4] , SNLS [5] and ESSENCE flBl assuming dynamical w 
with flat universe yields -1.38 < W0,n< -0-86 at 95% CL and w^^ = -1.12 ±0.13 at 68% CL. With additional BBN 
constraint of limit of expansion rate [IBill^l! ^0,0 = — E09±0.12 at 68% CL. The phantom field will finally dominate 
$H the universe in future, leading to Big Rip singularity [l8| . There have been many attempts to resolve the singularity 
5^ , from both phenomenological and fundamental inspirations . However fundamental physics of the phantom field is 
still incomplete due to severe UV instability of the field's quantum vacuum state [23| . 

This review interests in non-linear Schrodinger-type formulation of scalar field cosmology. We shall call the for- 
mulation, NLS formulation. In our NLS system, cosmological ingredients are scalar field and a barotropic fluid with 
constant equation of state, = w^p^. We also have non-zero spatial curvature. This is a system resembling of 
our present universe filled with scalar field dark energy and barotropic cold dark matter or of the early infiationary 
universe in presence of infiaton and other fields behaving barotropic-like considered in e.g. [11]. In such a model, the 
scale invariant spectrum in the cosmic microwave background was claimed to be generated not only from fluctuation 
of scalar fleld alone but rather from both scalar fleld and interaction between gravity to other gauge fields such as 
Dirac and gauge vector fields. 

Not long ago, mathematical alternatives to the standard Friedmann canonical scalar field cosmology with barotropic 
perfect fiuid, was proposed e.g. non- linear Ermakov-Pinney equation [2l|, |22|- Expressing standard cosmology with 
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fc > in Ermakov equation system yields a system similar to Bose-Einstein condensates [23|- Another example is 
a connection from a generalized Ermakov-Pinney equation with perturbative scheme to a generalized WKB method 
of comparison equation It was then realized that solutions of the generalized Ermakov-Pinney equation are 

correspondent to solutions of a non-linear Schrodinger-type equation, and then the NLS version of the Friedmann- 
Lemaitre- Robertson- Walker (FLRW) cosmology was formulated In the NLS framework, the system of FLRW 

cosmological equations: Friedmann equation, acceleration equation and fluid equation are written in a single non- 
linear Schrodinger-type equation. We will not prove it here but instead, referring to Ref . [25| . Few recent applications 
[26, 27, 28, 29] of the NLS formulation have been made and this review intends to conclude its major aspects. Note 
that application of Schrodinger-type equation to scalar field cosmology was previously made but in different form to 
tackle aspects of inflation and phantom field [s^l ■ 



II. SCALAR FIELD COSMOLOGY 



A. Friedman formulation 



We set up major concepts in this section before considering its application later. In the Friedmann system, barotropic 



fluid has pressure Pj and density with an equation of state, = w-yp-y = [{n~ 3)/3]p^ where n = 3(1 



Wj). Scalar 
= -1. The 



= w^Pff). To sum up, ptot = P-y + P4> and ptot — Pj+Ptp- Therefore n — means 
— —1/3; n — 3 for — 0; n — A for — 1/3; n — 6 for = 1. Barotropic fluid and scalar 



field pressure obeys p^ - 
others are: n = 2 for 

fluid are conserved separately. Dynamics of the barotropic is governed by fluid equation, pj = —nHp^ with solution, 
p,y = D /d^ , where a is scale factor. The dot denotes time derivative. H = a/a \s Hubble parameter and _D > is 
a proportional constant. Scalar field is minimally-coupled to gravity with Lagrangian density, C — (l/2)e0^ — ^{4') 
and is homogenously spread all over the universe. The scalar field density and pressure are 



P0 



1 



4>^ + ^(0) : 



1 



02 _ 



(1) 



The branch e = 1 is for non-phantom field case and e = — 1 is for phantom field case. Dynamics of the scalar field is 
controlled by conservation equation, e{(f> + 3H(j)) = —dV/d(f> , in which the spatial expansion H of the universe sources 
friction to dynamics of the fleld. The Hubble parameter is governed by Friedmann equation, = {k'^ /3)ptot — k/a^ , 
and by acceleration equation, a/ a = —{k? /Q){ptot + 3ptot) , which gives acceleration condition ptot < — Ptot/3. Here 



Ptot = WeffPtot, 

curvature and 



= SttG = 1/Mf,. G is Newton's gravitational constant. is reduced Planck mass, k is spatial 



Weff 



+ P-fW^f 
Ptot 



(2) 



If we express the field speed and the field potential in term of a{t) and time derivative of a(t), then 



e4>{tf 



H- 



nD 
3a" 



and = 



H 



2k 
3^ 



n — 6 \ D 
a" 



(3) 



B. NLS formulation 



NLS formulation is a mathematical alternative to the standard Friedmann formulation with hope that the new 
formulation might suggest some new mathematical tackling to problems in scalar field cosmology. In the NLS formu- 
lation, there is no such an analogous equation to Friedmann equation or fiuid equation. Instead both of them combine 
in single non-linear Schrodinger-type equation, 

u"{x) + [E- P{x)] u{x) = -!^u(a;)(^-")/" . (4) 
The links to cosmology are valid as one defines NLS quantities |25| , 

u{x)^a{tr^/\ E^-^D, P{x)^^a(tre4>{tf . (5) 
where ' denotes A/ Ax. Independent variable t is scaled to NLS independent variable a; as x = o'(t), such that 



x{t) = u{x) and (f>{t) = ip(x) , 



(6) 
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which gives e(j)(ty = eir ^'{x) . Hence e^'{xY = [A/ K^n) P{x) , and 



V'(x)=±^ [J^dx. (7) 



Inverse function i/j ^(x) exists for P{x) ^ and n 7^ 0. In this circumstance, x(t) = ip ^ o (j)(t) and the scalar field 
potential, V o a~^{x) and e0(i)^ can be expressed in NLS formulation as 



and 



The other quantities are 



12 , 2P n 12^; o 3fc 



(^.'f +44"' + ^"'/", (10) 



(U j H 2~2" 2" ' '11) 



Ptot = + (12) 

Ptot - -4^(«')' + ^"«"--^^i'^"- (13) 

= H = --uu", (14) 

71 71 

2Puu' + P'u^ ■ Uu'u f¥ , , 

-, ZB<i>= -V — • (15) 



K,\J Pen 

We shall see later examples that the program of NLS formulation must start from presuming the "wave function" , 
u{x) = a~"/2 = i{t)^ before proceeding to calculate the other quantities. We know that normalization condition for 
a wave function is |7/(a:;)pda: = 1 . If applying this to our NLS wave function, then x^Ax = 1 . In order to 
satisfy the condition, x must be constant (hence so is t) with an integrating constant = 1. In connecting Friedmann 
formulation to NLS formulation, we are forced to have u{x) = x{t). Therefore u{x) is, in general, non-normalizable. 



III. SLOW-ROLL CONDITIONS 



A. Slow-roll conditions: flat geometry and scalar field domination 

In flat universe with scalar field domination, H = —{k^ /2)<f>'^e . Hence for e — —1 (phantom field), 

0<aH^<'d, (16) 
i.e. the acceleration is greater than speed of expansion per Hubble radius, d/cH^^ and for e = 1 (non-phantom field), 

< a < aH^ . (17) 

Slow-roll condition [3l|, [l^l assumes negligible kinetic term, i.e. |e(/)2/2| <C ^(0) which makes an approximation 
~ K^V/S. This results in a condition \H\ <^ H^. Slow-roll parameter, e = ~H/H^ is hence defined from this 
relation. The condition \£<p/2\ ^ V{cj)) is then equivalent to |e| <C 1, i.e. — 1 ^ e < for phantom field case 
and < e <C 1 for non-phantom field case. Considering 77 ~ implying approximative constancy in H during the 
slow-rolling regime. For non-phantom field, this condition is necessary for inflation to happen (though not sufficient) 
[3^ 1 however, for phantom field case, the negative kinetic term always results in acceleration with < —1 then 
it does not need the slow-roll approximation. Another slow-roll parameter can be defined when the friction term 
dominates \4>\ ^ |3-ff0|. This gives the second parameter, rj = —(j)/H<j) and the approximation is made to I77I <C 1 
[3^ . The field fluid equation is then ~ —V^/ieH which implies that if e = — 1, the field can roll up the hill. With 
all assumptions imposed here, i.e. A: = 0, p-y = 0, \e(jP /2\ <C V and <C |3-ff(/)|, one can derive e = {l/2K^e) {V^/V)^ 
and 77 = {\/ K^){V(f,(f,/V) as known where the subscript denotes d/d0. 
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B. Slow-roll conditions: non-flat geometry and non-negligible barotropic density 



There are also inflationary models in presence of other field behaving barotropic-like apart from having only single 
scalar fluid [ssj . The the scale invariant spectrum in the cosmic microwave background was claimed to be generated 
not only from fluctuation of scalar fleld alone but rather from both scalar fleld and interaction between gravity to 
other gauge fields. Assuming this scenario with fc ^ and ^ 0, then 



K" ■„ k nn D 

The slow-roll condition becomes |K^e0^/6| ^ {n^V/i) — {k/a?) + (K^I?/3a") hence 

~ \ T VH\ 

3 3a2 18 a" 

implying | — {H/3) + (k/Sa^) — (nK^_D/18a")| <C H^- We can reexpress this slow-roll condition as 

\e + ek+eD\ < 1, 



(18) 



(19) 



(20) 



where = k/a?H^ and So 



^D/Qa^H'^ . Another slow-roll parameter 77 is defined as 77 = —(j)/H(j), i.e. the same 



as the flat scalar field dominated case since the condition ^ |3-ff0| is independent of k and p-y. 
Writing the condition |e</'^/2| < F in NLS form using Eqs. (P and 

j\ 2 



\P{x)\ « - 
n 



- +E 



4 



(21) 



If the absolute sign is not used, the condition is then e(p^/2 <C V, allowing fast-roll negative kinetic energy. Then Eq, 
pTjl . when combined with the NLS equation m, yields 

3u'2 



u" < 



n u 



kn 



(4 — Ti) /n 



(22) 



P' 




^ /6 -|- n\ u' 




< 






\ n J u 



Friedmann analog of this condition can be obtained simply by using Eq. ([3]) in the condition. Using Eq. (jlSp . the 
second slow-roll condition, \4)\ ^ |3i?(/i| in the NLS form is written as, 

(23) 



(24) 

(25) 
(26) 

(27) 



This condition yields the approximation 3He(f)^ ~ --dV/dif) which, in NLS form, is 

P' ^ 2u' 
~P ~ ~ 

The slow-roll parameters e, Sk and in NLS form, are 



2u' 



£k 



i^ku'^/" 
4u'2 



nE / M \ 2 



therefore 



£tot — £ + £k + £d — 



Hence the slow-roll condition, jetotl ^ 1, is just 



2 \u 



< 1. 



P{x). 



Another slow-roll parameter rj — —(j)/H(j) can be found as follow. First considering ^l^{x) = 4>{t) and Eq. ([7]), using 
relation d/dt ^ xd/dx, one can obtain 

n f u P' 



n f u lb" 
77 = - — — + 1 
' 2 \u' V/ 



The slow-roll condition \rj\ <C 1 in NLS form is just 



u P' 
it'2P 



2 \u'2P 



< 1. 



1 



(28) 



(29) 
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IV. ACCELERATION CONDITION 

For the phantom field, since its kinetic term is always negative and could take any large negative values, the slow-roll 
condition is not needed. The acceleration equation is taken as acceleration condition straightforwardly, i.e. a > 
hence 



(30) 



This, in NLS-type form, is equivalent to 



2 fu'V nk (u^l-" 
E-P > -— ' 



n \ u I 2 \ u 



(31) 



which is reduced to 



(32) 



with help of the Eq. (H]). Using Eqs. (I14p . the acceleration condition is just e < 1. 



POWER-LAW COSMOLOGY 



The power-law expansion a{t) = , with q > 1 is assumed here as the first step of calculation. In some high-energy 
physics models, during inflation, flat geometry and scalar field domination are assumed. The universe was driven 

by an exponential potential 1^((/)) — \q{?>q~ l)/{K^t1)\ exp | — K-\/2/g [(t>{t) — 0(to)]| [l^l- Also, at late time with 
dark matter component, the expansion could be power-law. Recent results from X-Ray gas of galaxy clusters put a 
constraint of q ~ 2.3 for fc = 0, g ~ 1.14 for k — —1 and q ~ 0.95 for fc = 1 [35|- For a flat universe, the power law 
expansion, a = t^, is attained when —1 < Wcs < —1/3 where q = 2/[3(l -I- Wcs)]. If using q = 2.3 as above, it gives 
= -0.71 (only flat case). Latest combined _WMAP5 results with SNI and BAO yield -0.0175 < ilk < 0.0085 
at 95% maximum likelihood [13] ■ The mean is flk = —0.0045 corresponding to closed universe with q = 0.986 [s^ . 
Assuming power-law expansion, the Schrodinger wave function is [1^ 



u{x) = x{t) = 



(33) 



Integrating the equation above so that the Schrodinger scale, x is related to cosmic time scale, t as 



X — x(t) — — hxo, and t{x) 

P 



1 



-P{x-x^f'f'' 



(34) 



where /3 = {qn — 2)/2 and xo is an integrating constant. The parameters x and t have the same dimension since (3 is 
only a number. Then the wave function is 



u{x) 



1 



qn + \ \ {x — xo) 



-1 qn/(qn-2) 



(35) 



which depends on only q and n. Wave functions for a range of barotropic fluids are presented in Fig. [TJ The result is 
confirmed by substituting Eq. (P5|) into Eq. The field speed and scalar potential are: 



2q 



2k 



nD 



kH^ 3*9" 



and V{t) 



g(3g-l) 



2k 



6\ D 



(36) 



From Eq. 



therefore the Schrodinger potential is found to be 



P{x) 



2qn 



{qn — 2)2 [x — xq)^ 



kn 
T 



-2 



{qn - 2){x - xq) 



2g(n-2)/(gn-2) ^2^2jj 



12 



(37) 



FIG. 1: u{x) versus x for power-law cosmology with q = 2. We set xo — 0. There is no real-value wave function for n = 3, 
n = 4 and n = 6 unless a; < 0. 




FIG. 2: P{x) plotted versus x for power-law expansion. We set q = 2,K = l,D = l and xo = 0. There is only a real-value P{x) 
for the cases k = ±1 with n = 4 because, when a; > 0, P{x) becomes imaginary in these cases. The physical value is when 
a; < since t has a reverse sign of x. 



With E ~ ~K n D/12, the Schrodinger kinetic energy is 



T{x) 



2qn 



1 



{qn — 2)2 (x — 



kn 

T 



-2 



{qn -2){x- xq) 



2q{n-2) / (qn-2) 



(38) 



A disadvantage of Eq. (|37|) is that we can not use it in the case of scalar field domination. Dropping D term in Eq. 
([57)1 can not be considered as scalar field domination case since the barotropic fluid coefficient n still appears in the 
other terms. The non-linear Schrodinger-type formulation is therefore suitable when there are both scalar field and 
a barotropic fluid together such as the situation when dark matter and scalar field dark energy live together in the 
late universe or in the inflationary models in presence of other fields behaving barotropic-like and single scalar fluid 
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[231 ■ P{x) is plotted versus x for power-law expansion with g = 2 in closed, flat and open universe in Fig. [2l One can 
check that the acceleration condition ([5^ for the power-law case is just g > 1. 

There is application of the NLS scalar field function i/j in Eq. ^ to solve for scalar field exact solutions in power-law, 
phantom expansion (a ~ (ta — O'j 9 < 0) ^-^d exponential (dc Sitter) expansion a ^ exp{t/T) f27l.l28|. For example in 
power-law case: 




2qn 



1 



'^V^ V ^(9"- ~ 2)^ (x — xof' 
The solution can be found only when assuming fc = 



kn 
27 



-2 



1 



[qn -2) [x ~ xo) 



2g(n-2)/(9n-2) 



12e 



dx . 



(39) 




2a 1- 



-f In 



2qn 



(40) 



When g = 2/n or n = 0, the field has infinite value, q and e must have the same sign for the solution to be real. The 
last logarithmic term does not restrict sign of q. This is unlike the solution obtained from Friedmann formulation 
which requires q < which violates power-law expansion condition {q > 1). Working in neither of them can obtain 
exact solution with fc ^ 0. In NLS formulation, we can not set D to zero while n is multiplied to the other terms 
then it can not be reduced to the scalar dominant case. This is a weak aspect. Obviously, the most difficult case is 
when fc ^ with D ^ 0. This case can not be integrated out in both frameworks unless assuming n — 2 (equivalent 
to w-Y = —1/3) which is not physical. 

There are other good aspects of the NLS formulation. Since transforming standard Friedmann formulation {t as 
independent variable) to NLS formulation {x as independent variable) makes n appear in all terms of the integrand 
and also changes fluid density term D from time-dependent term to a constant E, therefore the number of x (or 
equivalently t)-dependent terms is reduced by one and hence simplifying the integral ([7]). In the case of exponential 
(de Sitter) expansion using NLS formulation, the solution when k ^ and D can be obtained without assuming 
n value but n — 0,2, 3, 4, 6 must be given if working within Friedmann formulation. The phantom expansion case is 
very similar to the power-law case but only with different sign (see Ref. [11]). 



VI. PHANTOM COSMOLOGY AND BIG RIP SINGULARITY 



If we assume the expansion to a form, a{t) ~ (^a — t)'' with a finite time ia, one can see that q = 2/3(1 -I- Wcs) < 
(for a flat universe). This corresponds to Wes < — 1- Such equation of state is called phantom. The Schrodinger scale, 
X is related to cosmic time scale, t as 



x{t) 



[it. - t) 



and the wave function is 

u{x) = [f3{x - xo)]''"/'^"-' 
which is plotted in Fig. [3] with various types of barotropic fluid [2~ 



P(x) 



2qn 



{qn — 2)2 (x — xqY 



kn 
~2 



Therefore 

2q(n-2)/(gn-2) 



{qn ~ 2)(x - xo) 



12 



(41) 



(42) 



(43) 



Fig. m shows P{x) plots for three cases of fc with dust and radiation. P(x) goes to negative infinity at x = xq = 1. 
Expansion of the form, ait) ^ (ta — i)^ leads to unwanted future Big Rip singularity [ll|. The Big Rip conditions 
are that (a,ptot, |ptot|) ~^ oo which happen when t t~ in finite future time. Written in NLS language, if a — > oo, 
X ^ Xq and then it ^ 0+ (see Fig. [S]). Considering also Eqs. p2|) and p3|) . hence conditions of the Big Rip 
singularity are [1^ 



t - 




4^ 


X - 


~* ^0 


a - 


-> oo 




u{x) - 


-.0+ 


Ptot - 


-> oo 




u [x) - 


oo 


btoti - 


oo 




u'{x) - 


-> oo 



(44) 
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FIG. 3; Schrodinger wave function, u{x) when assuming phantom expansion. u{x) depends on only q, n and ta- Here we set 
~ 1.0 and q — —6.666. If fc = 0, g = —6.666 corresponds to = —1.1. 



Scalar field with dust (n = 3) Scalar field with radiation (/; = 4) 
P(x) P(_x) 





X 




X 












k = +1 






\ 





P(x) P(x) 




FIG. 4: Schrodinger potential in phantom expansion case for dust and radiation fluids with fc = 0, ±1. Numerical parameters 
are as in the u{x) plots (Fig. [3}. a;o is set to 1. For non-zero k, there is only one real branch of P{x). 



We have one less infinite value in NLS Big Rip condition, i.e. u{x) goes to zero. The NLS effective equation of 
state WcS = Ptot/Ptot can be expressed using Eqs. ([T^ and (fT5|) . Approaching the Big Rip, x ^ Xq ,u —>■ 0'^ , then 
WeS — > —1 + 2/3(7 where 9 < is a constant. This limit is the same as the effective phantom equation of state in the 
case k = 0. It is important to note that scalar field potential here is built phenomenologically based on expansion 
function, not on fundamental physics. 



VII. WKB APPROXIMATION 



WKB approximation in quantum mechanics is a tool to obtain wave function. However, in NLS formulation of 
scalar field cosmology, the wave function is first presumed before working out the shape of P{x). Procedure is opposite 
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to that of quantum mechanics. Hence the WKB approximation might not be needed at all for the NLS. Anyway, 
if one wants to test the WKB approximation in the NLS formulation, these below are some results. The WKB are 
valid when the coefficient of highest-order derivative term in the Schrodinger equation is small or when the potential 
is very slowly- varying. Consider linear case of Eq. {k = 0), 



[P{x) - E]u = 0. 



(45) 



In Figs [2] and m the left-hand side of P{x) is physical since it corresponds to positive time. In most regions, there are 
ranges of slowly varying P{x) at large value of in which the WKB is valid. The approximation gives 



a ^ A exp 



± 



{2/n)i / ^/E^^¥(x)da 



(46) 



where A is a constant. 



VIII. CONCLUSIONS 



Here we conclude aspects of NLS-type formulation of scalar field cosmology. The NLS-type formulation is well- 
applicable in presence of barotropic fluid and a canonical scalar field. There are few advantages of the NLS formulation 
as well as disadvantages to the conventional Friedmann formulation. With hope that some more interesting and useful 
features could be revealed in future. 
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